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NEW PHASE TRANSITIONS IN CHERN-SIMONS MATTER THEORY 


ALI ZAHABI 


ABSTRACT. Applying the machinery of random matrix theory and Toeplitz determinants we 
study the level k, U ( N ) Chern-Simons theory coupled with fundamental matter on S 2 x S' 1 at 
finite temperature T. This theory admits a discrete matrix integral representation, i.e. a unitary 
discrete matrix model of two-dimensional Yang-Mills theory. In this study, the effective partition 
function and phase structure of the Chern-Simons matter theory, in a special case with an ef¬ 
fective potential namely the Gross-Witten-Wadia potential, are investigated. We obtain an exact 
expression for the partition function of the Chern-Simons matter theory as a function of k, N, T, 
for finite values and in the asymptotic regime. In the Gross-Witten-Wadia case, we show that 
ratio of the Chem-Simons matter partition function and the continuous two-dimensional Yang- 
Mills partition function, in the asymptotic regime, is the Tracy-Widorn distribution. Conse¬ 
quently, using the explicit results for free energy of the theory, new second-order and third- 
order phase transitions are observed. Depending on the phase, in the asymptotic regime, Chern- 
Simons matter theory is represented either by a continuous or discrete two-dimensional Yang- 
Mills theory, separated by a third-order domain wall. 


1. Introduction 

Study of the phase structure of (p + 1)— dimensional Yang-Mills (YM) theory coupled to 
matter on S p sphere and especially the confinement/deconfinement transition in large N gauge 
theories on sphere sheds light on various phenomena in gravitational phase transition such as 
black hole nucleation, via the AdS/CFT correspondence. In this direction, the phase structure 
of Chern-Simons (CS) theory coupled to matter and its gravitational dual, Vasiliev higher spin 
gravity, flVal . have been studied recently, UCh-Mil and llGUl . In a recent work fTall . it has been 
argued that level k, U(N) CS theory coupled to fundamental matter on S 2 X .S’ 1 at finite tem¬ 
perature T, can be explained by a discrete unitary matrix model, and its phase structure, the 
eigenvalue density in different phases, has been studied in the limit of large k. N. T, with fixed 
parameters A = N/k. £ = T/N, by using steepest descent method. This study is based on 
the reduction of the Chern-Simons matter (CSM) theory first to the CS theory with an effective 
potential and second to a discrete version of two-dimensional YM theory and its representa¬ 
tion, one matrix model. In a special case with Gross-Witten-Wadia (GWW) effective potential, 
IIGr-Wift and IWal . in addition to well-known lower-gap phase, an upper-gap phase is observed, 
flU and IITakl . in which the upper bound of the eigenvalue density is saturated. 

In another discipline of research, discrete Toeplitz determinants are studied in itBa-Lil . and 
by using the continuous orthogonal polynomials their precise forms are obtained in terms of the 
continuous Toeplitz determinants and the Fredholm determinant. A key fact to understand the 
application of the Toeplitz determinant in the CSM theory is the Heine-Sezgo identity which 
relates the matrix model integral and the Toeplitz determinant, see appendix A. 

Using the relation between Toeplitz determinants and matrix integrals, and the above results 
for discrete matrix model representation of the CSM theory, in this paper, the techniques and 
methods from discrete Toeplitz determinant is employed to study the CSM theory. The matrix 
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integral representations of the YM theory and CSM theory provide the application of well- 
developed techniques of random matrix theory and specially the Toeplitz determinant. 

In our setting, the partition functions of YM theory and CSM theory with an arbitrary poten¬ 
tial V(T) can be written as Toeplitz determinant and discrete Toeplitz determinant with respect 
to the probabilistic weight f(T) = e ' <T K respectively, 

(1.1) Z YM :=Z N (V) = B N (f ), Z CSM -.= Z k N {V) = V k N {f). 

In this study, we move towards an exact determination of the thermal partition function of 
CSM theory. In fact, by adopting recent results, HBa-Lil . about the discrete and continuous 
Toeplitz determinants, a new explicit expression for the partition function of CSM theory in 
terms of the partition function of YM theory, as a function of k, N, T at the finite values and 
in the asymptotic regime, is obtained. Especially, a CSM theory with GWW effective potential 
is studied in detail and analytic forms of the effective partition function and free energy of 
CSM theory are obtained. Furthermore, the sub-leading corrections in the asymptotic results 
are computed and the result is interpreted as the energy of the upper bound. These results lead 
to a phase transition in the asymptotic regime of the discrete matrix model, related to the phase 
transition between upper-bound phase and no-upper bound phase. 

Although, the effective GWW potential is not a microscopic potential and it does not exactly 
characterize the CSM theory, but the CSM model with GWW potential shares some features 
with the real CS theory such as the minimum of any effective potential v(U) (as a function of 
holonomy U) is always at U = I and the depth of the potential increases as a function of (. 

In brief, the phase structure of the model in the asymptotic limit of large N, T , k is studied via 
the asymptotic behavior of the ratio of partition functions which is obtained as the Tracy-Widom 
distribution. Then, from analytic expression for the free energy, obtained from the expansion 
of Tracy-Widom distribution, new second-order and third-order phase transitions in this model 
are determined. Moreover, relation between the obtained result in this paper and other recent 
results is explored. 

The phases of CSM theory, can be discussed from the viewpoint of the existence of the con¬ 
tinuum limit. In fact, the continuum limit does not exist in whole moduli space of the parameters 
and in fact, N, T, k should satisfy certain relations in order to determine the continuum limit. 

This paper is organized as follows: In section two, the CSM theory and the discrete matrix 
integral representation of the partition function are explained in brief. In section three, the main 
results of this study are expressed. First, the partition function and free energy of the CSM 
theory with GWW potential are explicitly obtained and subsequently the phase structure of the 
model with GWW potential is studied. In section four, possible interpretations of the obtained 
results are discussed and directions for further studies are introduced. Finally, there are two 
appendices, where necessary materials from the random matrix theory such as orthogonal poly¬ 
nomials, Toeplitz determinant and Riemann-Hilbert problem (RHP) and their inter-relations 
are introduced and well-known results in two-dimensional YM theory with GWW potential in 
terms of Toeplitz determinant and RHP problem, and some results about the partition function 
and phase structure in the CSM theory are summarized. Finally, some necessary definitions and 
results about Tracy-Widom distribution are collected. 

2. Chern-Simons matter theory 

In this section, the level k, U(N) CSM theory coupled to fundamental matter at temperature 
T on S 2 x .S' 1 is described. Especially, the CS theory coupled with the matter in fundamental 
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representation on S 2 x S 1 , in T, N, k —> oo limit with the fixed’t Hooft parameter, A = N/k, 
and ( = T/N, will be studied. In this part, we explain two relevant and necessary results from 
the literature; I) the action of pure CS theory coupled with fundamental matter in the’t Hooft 
limit and at high temperature, up to leading order in N, can be expressed as the sum of the action 
of pure CS theory and an effective action, i.e. effective potential which is a local function of 
holonomy matrix U(x ) with eigenvalues e ia around S 1 with x £ S 2 . II) The partition function 
of the CSM theory with effective potential can be written as a discrete version of the matrix 
integral representation of the 2d YM theory with the effective potential. These observations are 
explained in this section and in the Appendix A, not only for the article to be self contained but 
also because of then - importance and usefulness in understanding the physics and mathematics 
of the problem studied in this article, since our new results are based on these observations. 


2.1. CSM theory and discrete YM theory. CS theory is a three-dimensional topological field 
theory of the Schwarz type ffScl . with no metric dependence in the action. Consider a three- 
dimensional manifold M with a principal G —bundle with a connection A, then the CS action 
is defined by 

(2.1) Scs(A) = ^-[ Tr(A A cL4 + A A A A), 

4vr J M o 

where k is called level and it is an integer in the quantized theory. The above action functional 
defines a three-dimensional topological field theory, i.e. a field theory with observables, explic¬ 
itly independent of the metric on M. The partition function of the theory is given by a path 
integral, 


( 2 . 2 ) 


Zcs(A) = / [DA] e 


i Scs 


and the observables of the theory are defined via the correlation functions of the Wilson loops, 

WM 1 


W Ci ,rM) = Tr 


in Ri representation with closed loops Ci, as, 


(2.3) 


ilW Ci ,RM) ) = / [DA] e iScs [UWc^rM) 


\i=l 


\i=l 


In this study, we consider G = U(N). The partition function and correlation functions of 
the theory define global topological invariants of the manifold as well as other topological 
invariants associated to the manifold such as knot invariants, Jones polynomials, etc. llWill . 
In this study, we are interested in CS theory on curved background S 2 X S 1 . Such curved 
backgrounds provide enough symmetry for explicit computation of the partition function and 
correlation functions. For example, the partition function of CS theory on Seifert manifolds 
such as S 3 admits a matrix model representation and thus analytic closed expressions for the 
partition function have been obtained in llMalll . For a general review on this and other different 
aspects of CS theory, see llMa2l . 

The gauge theory in this study is a level k, U(N) CS theory at finite temperature and finite 
N, as well as its high temperature and large N limit, coupled to fundamental matter fields. 
The behavior of the CSM theories and then - partition functions and free energies as functions 
of £ and A has been a subject of intense studies during past years. In the limit A —> 0, the 
behavior of CSM theory and its thermal free energy has been studied in USh-Yil . ifGTl . UCh-Miil . 
Fundamental matter CSM theories at finite ’t Hooft coupling A and at large N have been studied 

3 












in m, fiCh-Mil . lAh-Gu-Yali . HAh-Gu-Ya21 . llMa-Zhll . HMa-Zh21 . IIBa-Hel . Hi, (SI, ffiCiul . 
flAhll . The phase structure and thermal free energy of some examples of such theories have been 
studied in IfTal . I'Takl . 

The partition function of the CSM theories contain the CS action and a matter action which 
shows the interaction with the matter fields, 

(2.4) Z CS m(A,h) = I [VA\[Vfi\ e*cs-s^ 

where p represents the matter fields and <S /( is the coiTesponding action. The matter content 
of the CSM theory consists of vector massive fields in the fundamental representation such as 
fermionic and bosonic fields minimally coupled to CS fields, for further details see references in 
the previous paragraph. In this study we will not discuss any computations containing explicit 
forms of the matter fields, thus we will keep to our general form of matter action and we will 
not be more explicit about it. As we will see in the following, in the high temperature limit, the 
matter fields have masses of order the temperature and the effective action only depends on the 
holonomy field. 

The goal of this part is to compute the partition function of CSM theory. In order to do 
this computation, roughly speaking, we have to sum all the vacuum graphs of the massless and 
massive fields and thus the partition function of CSM theory is obtained by summing up all 
the separate contributions from the vacuum graphs of massive and massless fields including the 
holonomy field. In brief, we first evaluate the path integral on the matter fields which gives us 
the partition function as a path integral over massless fields with a CS action and an effective 
action for the holonomy field. The effective action is directly obtained from path integration 
over the massive fields. More precisely, after integration over the massive fields, a local effective 
action, S e ff(U), for the 2d unitary matrix valued holonomy U(x) around 5 1 with x £ S 2 is 
obtained and the CSM partition function takes the following form, 

(2.5) Zcsm(A) = j [DA] jScs-S eff ( u ) 

Here, we briefly discuss the behavior of the effective action such as its series expansion and 
its leading term in the large N limit which is a local function of holonomy field. The key fact 
in evaluation of the effective action and CSM partition function is the following observation in 
flAhl . In order to obtain the effective action, we integrated out the vacuum graphs of the matter 
fields and since these matter fields develop large thermal masses, thus S e f f(U) can be expanded 
as the following local series, 

(2.6) S eff = j d 2 x{Tv{U{x )) + (U)TrDiUD i U + ..., 

where v(U(x)) and v\ (U(x)) are local effective potentials which their actual forms depend 
on the matter content and CS couplings. In the large N limit, the explicit computation of the 
effective potentials as a function of A, for any given vector matter CS theory is feasible. For 
example, for CSM theory with U = I, and minimally coupled fundamental fermions to the CS 
fields, the effective potential computed by using large N techniques in the light-cone gauge, 
(Gil . This computation has been generalized to other theories and different holonomy matrices, 
iDail . liAhl . IITall and references therein. It has been argued in itAhl . that the effective potentials 
are of order N, and since T = (N, the first term in the expansion of the effective potential is 
of order N 2 and the second one is of order N and the rest are further suppressed at large N. 
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Therefore, at large N with T ~ N, up to leading order in N, the entire effect of matter loops 
in vacuum graphs on the dynamics of the CS theory is represented by the high temperature 
effective action of the following form, 

(2.7) Seff = T J d 2 Xy/gv(U(x)). 

And finally the thermal partition function of the CSM theory becomes 

(2.8) ZcSM = J[VA]j&fMTr(AdA+%A*)-Tl'P xV gvmx)) =< e ~TVMU) ^ 


where V 2 is the 2 d volume of S 2 and since the pure CS is topological, the expectation value of 
the exponential of the integral of the effective potential over S 2 is independent of x. Therefore, 
the CSM partition function can be written as a topological observable of the pure CS theory, 
i.e. an expected value of linear combination of the Wilson loops in pure CS theory. For more 
convenience assume V 2 = 1. 

An explicit, matrix integral representation of the CSM partition function (2.8) can be ob¬ 
tained by using the path integral techniques, iThl . | |B1-Thl l. see Appendix A. Roughly speaking, 
after gauge fixing and evaluating the integral of the CS action over S 2 , the path integral over the 
space of gauge fields reduces to a matrix integral over the eigenvalues of the holonomy matrix. 
This matrix integral representation of Zqsm is a discrete version of the matrix integral represen¬ 
tation of the two-dimensional YM theory partition function, Eq. (A. 1 1 with Vym = NQu(U). 

is given by a 


As a result of this computation, partition function of CSM theory, Eq. 
discrete unitary matrix integral and equivalently, by using Fleine-Szego identity, as a discrete 
Toeplitz determinant, 


ZcSM = j [VA\ e iScs-T f d*xygv(U(x)) 


(2.9) 


N 00 


n E [II ( 2 sin( a,( ^ } ~ ^M V^)] 


i= 1 rij=—00 la 


= det 


l E 


z£d\ 


N -1 


3 , 1 =0 


= D k M(h), 


where rq’s are integer and a, (Tt) = are the eigenvalues of the holonomy matrix, (dis¬ 
tributed on the unit circle with distance 2n/k between two consecutive eigenvalues), and z = 
e ia and d\ = {z e <C\z k = 1} is a finite discrete subset of a unit circle S 1 with size |g?i| = k 
and / = e~ N< * v ( u } is the weight function. 

A toy model of the CSM theory, Eq. ( |2.9[ ), with interesting features is the CS theory with the 
GWW potential, v(U) = — \Tr(U + W), or v(a) = — £Y =1 cosa*. In the following section, 
based on the above Toeplitz determinant formulas, exact analytic formulas for the free energy 
of the CSM theory with GWW potential are obtained. Furthermore, using the obtained explicit 
results for the free energy, the phase transitions and their orders between different phases of the 
theory are studied. 
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3. New results in Chern-Simons matter theory 

In this part, the mathematical results in Toeplitz determinants have been adopted and trans¬ 
lated in the language of gauge theory and their deep implications and applications in the study of 
gauge theory, e.g. in determining the free energy of the theory and revealing the phase structure 
have been studied. As we will see, this new method in study of gauge theory is not only able to 
reproduce a very recent important result about the phase structure of CSM, produced by other 
plausible methods in IITall and ttTaktt . but also provides explicit expressions for the free energy 
of the gauge theory which contains important new results for the order of phase transitions in 
gauge theory. 

In this section, a new explicit results for the free energy and new phase transitions in the CSM 
theory with GWW potential are obtained. In order to obtain the final result, in the first part, the 
partition function of the CSM theory with an arbitrary potential is obtained via a recent result 
in the context of the discrete Toeplitz determinants. In the second part, the CSM theory with 
GWW potential is studied and explicit asymptotic results, in the limit of large parameters, for 
the partition function is obtained. In the third part, a careful study of the obtained results leads 
to an explicit expression for the free energy of the CSM theory with GWW potential which 
consequently reveals new phase transitions in this case. In the fourth paid, some consistency 
checks of the obtained results as well as some explicit computations for the free energy in 
different limits are performed. 


3.1. CSM theory partition function. It has been recently shown that the analysis of discrete 
Toeplitz determinants can be done by means of continuous orthogonal polynomials associated 
to continuous weight function, HBa-Lil . The main result in that study is the explicit relation 
between the discrete and continuous versions of the Toeplitz determinant for a given weight 
function in terms of a Fredholm determinant. In this study, we investigate the implications of 
those results for gauge theory. In other words, via the relation between the Toeplitz determi¬ 
nants and partition functions, see appendix A, the relation between the discrete and continuous 
Toeplitz determinants is used to connect the partition functions of CSM theory (discrete YM 
theory) and continuous YM theory. The physical meanings and implications of this connection 
for phase structure of the CSM theory are studied. 

In order to make the connection between Toeplitz determinants and CSM gauge theory pre¬ 
cise, we studied the relations between the parameters in both subjects. In both subjects, there 
are parameters for indicating the features of the systems and these parameters are closely related 
by their meanings in each subject. For example, the rank of the gauge group in CSM theory is 
equivalent to the size of the Toeplitz determinant and the level of the CSM theory is equivalent 
to a parameter in discrete Toeplitz determinant that shows the discreteness. The temperature of 
CSM theory is an extra parameter in the weight function of the Toeplitz determinant. In sum¬ 
mary, there is a one to one correspondence between the parameters in CSM theory and Toeplitz 
determinants which can be easily observed by comparing two subject. 

After adopting the results in discrete Toeplitz determinants for the gauge theory in this study, 
the fol lowing explicit formula for the partition function is obtained from Eq. ( A.9[ > and Eq. 
(A. 10) in Appendix A. In our setting, z = e 1Q and domain d\ = {z € C| z k = 1} is a finite 
discrete subset of a unit circle S 1 with size \d\\ = k (level of CS theory plays the ro le of 
discreteness in Toeplitz determinant) and f(z) = e v ' z ' 1 is a weight function, then, Eq. (A.9i 
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and Eq. (A.10) imply 


(3.1) 


Z-CSM 

Zym 


D^(/,di) 


D N (f) 

where K is the integral operator with a kernel. 


= det (1 + A'), 


,,.\ -nPn(z)p* n (w) -p* N (z)p N (w) 

(3.2) K(z,w) = z - - \ Vk{z)v k (w)e v ( z >e v ( w >, 

1 — z 1 w v 


and pn{z) is an orthogonal polynomial with respect to weight function e v ( z \ (see Appendix 
(A.2)) and p* N (z ) := z A p J v(^ _1 ) (on the circle S 1 , z _1 = z), and by using 7 (z) = z k — 1 in 
our domain, discrete function v k {z) is obtained from Eq. (A. 12 1 as, 

Vk(z) = 


(3-3) 


1 -z k 
~ — k 


l-z- 


z € -Si 

-2 e -SL 


Finally, the Fredholm determinant for the integral operator K is formally defined by, 

' (-l) r 


(3.4) 


det(l + K) = 1 + 


n= 1 


n\ 


, J det (K(zi, Zj))ij = 1 dzi...dz n . 


Up to a normalization, the free energy of CSM theory is obtained by taking the logarithm of 
Eq. ((3 3 , 

(3.5) Tcsm = ?ym + log det (1 + K). 

For a trace class operator K, we have, 

TrK n 


(3-6) 


Fcsm = Fym + 2>(log (1 + K)) = Tym + > v ( _ l) 


OO 

D 

n =1 


\n+1 • 




The above formal expression ( |3.5[ ) is a new exact result for the free energy of the CSM theory at 
any values of N , k, T. In principle, it contains more information about CSM theory in compari¬ 
son to the existing results in the literature which only discuss the free energy in the large N,k,T 
limit. In fact, its explicit evaluation reduces to evaluation of the Fredholm determinant with a 
kernel associated to the effective potential of the CSM theory. Moreover, Eq. ( |3.5[ ) basically 
contains explicit information about the phase structure of the CSM with any potential described 
in Eq. (2.9 1 . Roughly speaking, the kernel K which depends on the weight function / contains 
information about the free energy and it determines the phase structure of the theory. 

In fact, as it will be clear in the following section in the example of CSM theory with GWW 
potential, by using this new method we can explicitly obtain the free energy and the phase 
structure of the model. And since the CSM model with GWW potential shares some features 
with the real CS theory we expect that CSM theory with GWW potential has the main qualitative 
features of the phase structure of realistic CSM theories with other different potentials. The 
CSM theory with any potential and its phase structure can be exactly studied by using the 
explicit form of the free energy given by Eq. (3.5). However, in this study we focus on the 
toy model with GWW potential as the first step towards CSM theory with other complicated 
potentials. 

In gauge theories, we are interested in the limit of large parameters, in which the ratio of 
the partition functions can be studied via the asymptotic analysis of the Fredholm determinant. 
This analysis is mainly based on the asymptotic analysis of the orthogonal polynomials, 
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IIMa-Mcfl . In this study, instead of repeating the standard asymptotic analysis of the Fredholm 
determinant we rather prefer to focus on the physical aspects of the obtained results from these 
analysis. As the simplest example of such analysis, we consider a case of a positive constant 
weight function, f c = const. In this case, the asymptotics of orthogonal polynomial as N 


oo, 


(3.7) 


, is 


PN (z) = 


z N 0(e ~ CN ) 
„—CN\ 


<D(e~ CN ) 


\z\ > 1 + e 
\z\ < 1 - e ’ 


where C is a constant, and by using Eq. (3.31, the discrete function Vk{z) for large k becomes 
(3.8) 


„*(.)< 2 M\ 

V ' “ | 2(1 + e)~ k 


\z\ = 1 -e 
\z\ = l + e 


oo ) for the ratio of partition functions, Eq. (3.1), 


Above asymptotic formulas for orthogonal polynomials and discrete function, imply the fol 
lowing asymptotic results (in the limit N, k 
with constant potential, 

r(fc,di) 


a'k 


= 1 + 0{e~< k+N) ), 


(3-9) lim 

k—N—>oo,N —>-oo 2'N\Jc) 

where c is a positive constant and we have used Eq. <o>- The analysis details of the proof of 
the above result can be found in I Ba-Lil . We are interested in the physical meaning of the above 
result. In fact, Eq. ( |3.9[ ) is expected, because, when the discrete structure vanishes in the limit 
k —> oo, then, in the leading order, the partition function of the continuum YM theory is the 
continuum limit of partition function of the discrete YM theory. In other words, for a constant 
potential, we expect free energy be a smooth function in the asymptotic limit. The continuum 
limit is smooth in the sense that we do not expect any discontinuity in the behavior of CSM 
theory that leads to a phase transition in the continuum limit. However, this soft continuum 
limit is not always the case and as it will be studied in this paper, the asymptotic limit of the 
ratio ( |3.1[ ) depends on the weight function in a nontrivial way and the continuum limit is only 
obtained in a special region of the moduli space of the parameters. 

3.2. Partition function and phase structure of CSM theory with GWW potential. In this 
part, we explicitly demonstrate the usefulness of Eq. ( |3.1[ ) in gauge theory and its ability to 
produce physical results in CSM theory. In an important example, the partition function of 
CSM theory with GWW potential is obtained explicitly in terms of the partition function of 
the GWW model, (see Appendix A. 1.2), and by careful analysis of the asymptotic regime, new 
phase transitions arc traced. 


In the finite regime, the following equation, obtained from Eq. (3.11, determines the CSM 
partition function with GWW potential, 

(3-10) 4sm W) = 4°I W) det (1 + K GWW (z, w)), 

where Kqww(z,w) in principle can be obtained from Eq. (|3.2|) by putting Vaww( z ) = 
— ^(z + z -1 ) and using appropriate pn(z) obtained from Eq. (A.261 with Vcww , see ifBa-Lil . 
and v(z) from Eq. (3.31. 

In the limit N, T, k —> oo, asymptotic results for the ratio, Eq. (13. 1 1), in the case of GWW 


potential can be obtained either from a direct asymptotic analysis of the Fredholm determinant, 
similar to the analysis for the case of the constant potential in the previous part, or by using the 















adopted version of the results for Toeplitz determinants in the context of Brownian motion, see 
Appendix A. 1.1. In this study, we will not repeat any details of such analysis in the Brownian 
motion, instead we formally use the results in Brownian motion as some mathematical results 
and we further explore their implications for gauge theory. As a mathem atical result, in domain 
d s = {z € C| z k = s}, for fcww = e~ Vaww by putting together Eq. (A.131 and Eq. (A.14), 
we obtain, 

I Z%Ugww, d s ) 

f\s 1=1 ZnUgww) 27ris 


(3-11) 


lim 

N,T,k—*oo 


= F(^), 


a 


where F is the Tracy-Widom distribution (see appendix B) and functions // and a are defined 
by 


(3.12) 


H : = 


S N 

(2 y/NT N 


+ T N >T 


< T 


i „ i 
3 T 3 


a : = 


3 'T 1 3 ( A — -I- A / — 

61 J u/ 'r \i NJ 


N >T 
N <T 


For s = 1, Eq. ( |3. 1 1 [ ). implies the following new and useful result in gauge theory which is 
an explicit relation between the partition functions of CSM theory and YM theory with GWW 
potential, 

r(GWW) 


(3.13) 


lim 


■'CSM 


N,T,k ^00 7 (GWW) 
^YM 


= n 


k — /r, 


a 


As a remark, by using an expression for the Tracy-Widom distribution in terms of the kernel of 
the Airy function, F = det(l + Kait), (see Appendix B), Eqs. (3.10) and (3.131, imply that 

li m iV,T,fc-»oo Kgww(z, w) = K Ai (z, w). 

In the following, we extract the physical meaning behind the simplicity and elegancy of the 
above mathematical result, Eq. 


of the partition functions, Eq. 


S . In the next step, an asymptotic formula for the ratio 
is obtained via the asymptotic formula for the Tracy- 


Widom distribution. Let us denote x = then it is easy to see that the asymptotic limit 
k,N,T H> 00, leads to x -X ±00 for k > fi and k < fi, respectively. The asymptotic formula 
for the Tracy-Widom distribution, Eq. ( B.5| ) in Appendix B, is 

_-.x3/A 


(3.14) 


1 - 0{e 
0 (e-N 3 ) 


F(x) = 

Let us denote the ratio by lZ(N,T,k) = linrv,T.fc- 


x 

x 


00 
—00 

7 (GWW) 
-'CSM 


■OO 7 (GWW) 


= lim x ^ ±oc F(x), then by 


using Eq. (3.14 1 and putting T = (N, k = A 1 At, the ratio of the partition functions is given 

by. 


(3.15) 


K(N, T, k) = 


1 -0(e~ N ) 
0 + 0(e~ N2 ] 


k > (1 + e)fi(N, T) 
k <{l- e)n(N,T) ’ 


where e is a positive infinitesimal parameter. Equivalently, the above result, up to leading order, 
can be written separately for two phases, 

(3.16) 


for N > T : 1Z(k, N, T) 


k > N + T 
k<N + T 


forN < T : H{k,N,T) 


k > 2 y/NT 
k < 2 y/NT 






















Remembering that the n— th order of phase transition is defined by a discontinuity in the 
n —th derivative of free energy, thus, since the ratio of the partition functions jumps across 
the line k = p, one would expect a phase transition of order zero, however, more careful 
analysis by considering the sub-leading terms in the ratio is needed to determine the order of 
the phase transition via the explicit form of the free energy. As it is observed, CSM theory 
has a complicated moduli space of three parameters N, T, k indicating a phase transition at 
asymptotic regime. 

In order to explain the above result for the ratio, we need to understand the phase structure 
of the eigenvalues of the holonomy matrix, see Appendix (A. 1.2). By definition, density of 
the eigenvalues of the holonomy matrix is positive and normalized by p(a)da = 1. The 
eigenvalue density is obviously restricted from below by zero lower bound. The lower gaps 
in the eigenvalue density are defined by {a\a ~ a + 2tt,p(o) = 0} and this leads to the 
lower-gap phase of the CSM theory. As it is explained in Appendix (A. 1.2), and in Ail , in 
the CSM theory the eigenvalues of the holonomy matrix are discrete and distance between two 
consecutive eigenvalues is 2-jr/k. The discretization of the eigenvalues leads to an upper bound, 
^ X = 2 ^t, in the eigenvalue density, as 


(3.17) 


0 < p(a) < 


2vrA 


The upper bound on the eigenvalue density leads to the upper-gap phase with upper gaps de fined 
by {a\a ~ a + 27t, p(a) = The phase structure of the CSM theory defined by Eq. (2.9), 


can be understood via the comparison of the competing forces in the theory. There are two 
competing factors in this problem: I) the attractive potential V(U), of order 0(N 2 ) at leading 
order with T = (N, which tends to clump the eigenvalues aij’s and II) the repulsive force, of 
order 0(N 2 ), from the measure VU or equivalently from the Vandermonde determinant which 
repels the eigenvalues, roughly because it vanishes when two eigenvalues coincide. Roughly 
speaking, competition of the attractive force from the potential and a repulsive force from the 
Vandermonde determinant in the matrix integral, results to a phase transition. In fact, there 
are two phases; in one phase the attraction force is weaker than the repulsive force and the 
eigenvalues has support on the whole circle. In another phase, the repulsive force is weaker and 
this leads to the existence of a gap in the eigenvalue density and therefore the eigenvalues have 
only support on finite arcs on the circle. 

In summary, the saddle point analysis in the variational problem ( |A.16 ) indicates that in the 
CSM theory with the GWW potential, there are four phases, namely no-gap phase, lower-gap 
phase, upper-gap phase and two (lower and upper)-gap phase. Depending on the parameters A 
and C the phases of the CSM theory are classified in IITal . For summary of details of the phase 
structure of this model see the Appendix (A. 1.2) 


Finally, a possible interpretation of the above result (3.16) can be expressed in terms of the 


continuum limit of the CSM theory and/or existence of the upper-gap phase. Notice that the 
difference between the CSM theory and the YM theory is the presence of the discreteness and/or 
upper bound in eigenvalue density for the CSM theory. Therefore, the continuum limit of the 
CSM is expected in the regime or the phase that upper bound is not saturated. On the other 
hand, in the continuum limit the partition function of the CSM theory becomes the partition 
function of the YM theory and the ratio becomes one. Therefore, as we are interested in the 


phase structure of the theory at large N. k. T with fixed (, A, from Eq. (3.16) with the change 
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of variable in terms of A and £ the ratio with the correct interpretation can be written as 


(3.18) forC<l: n(X,0 


1 no upper-bound/continuum phase , A 1 > £ + 1 
0 upper-bound/discrete phase, A" 1 < +1 


(3.19) for(> 1 : 


1 no upper-bound/continuum phase, A _1 >2v^ 
0 upper-bound/discrete phase, A" 1 < 2- v /( : 


In the next part, explicit calculations for the free energy in each region of the moduli space 
determine the order of this new phase transition and other features of the CSM theory. 

Above interpretation is consistent with the summarized results in the Appendix (A. 1.2). This 
consistency can be seen by comparing our above result with the rearranged form of (|A.24|), 
(3.20) 


forX< \ 


no-gap or lower-gap, 
two-gap, 


£ < 4A 2 

£ > IX 1 


for A > ^ : 


j no-gap, C < 

1 upper-gap or two-gap, Q > 


However, in contrast to Ifm . our phase structure result is completely based on the analysis of 
the ratio of partition functions, without calculating the eigenvalue density. Moreover, as it will 
be clarified, our result contains more information about the phase structure such as the explicit 
formulas for the free energy which gives the orders of phase transitions between continuous and 
discrete YM theories (CSM theory) and also between discrete YM theories (CSM theories) in 
different phases. The actual properties and features of these new phase transitions can only be 
described by careful analysis of the Tracy-Widom distribution which will be presented in the 
next part. 


3.3. Free energy of CSM theory and orders of phase transitions. In this part, an explicit 
formula for the free energy of the CSM theory in different phases is obtained. The free energy 
determines the order of the new phase transition introduced in the previous part. The obtained 
results will be checked in nontrivial calculations. Furthermore, some limits of the free energy 
in interesting points of the moduli space are computed. 

In order to study the phase structure of the theory, the first step is to compute the sub-leading 
corrections to ratio formula, Eq. (3.151, by using the expansion of the Tracy-Widom distribu¬ 
tion. In the asymptotic regime, the free energy of the YM and CSM theories is defined via their 
partition functions, Tym/CSM = lirn.v-.oc 7 ^ log Zym/CSM- I n the case of GWW poten¬ 
tial, by using Eq. (3.131, and the integral representation of the Tracy-Widom distribution (see 
Appendix B), the free energy of the CSM theory is given by 


1 1 

(3.21) Fcsm = Fym + = Fym - ^2 [J (s - x)q 2 (s)ds\, 

where x = In different phases, by using YM free energy, || Gr-Witt . 


(3.22) 


-r-(GWW ) 
■'"YM 


(0 
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0 < c < 1 
C> 1 











the free energy of the CSM theory is given by 
(3.23) 

FCSM = TYM + ^ l ° gF(N l) = | C - I - ^ - NhXHsM C > 1 ' 


where 


(3.24) 


] ■= 


k—(N+T) N -% _ A"!—(C+l) 

l l ■‘ v •’ — I I 

2"3T3 2“3£3 


fc—(2y r /VT) 


2-§T3(t/fW$)S 


T N~ 3 = 


A -1 —2£2 
“ 2 —T - T-T 7 T 


2-3C3(C2+C“2)3 


C< 1 
C> 1 


Following the interpretation of the ratio result in the previous section, Eq. (3.21 1 can be written 
as 


(3.25) 


Fcsm = FyM + Fub, 


where T u \, = log F(N3 j) = — J'^% (s — N3 j)q 2 (s)ds is the upper bound energy. For 

j < 0, in the large N limit, by using the expansion of q(s) in the limit s —> ±oo, (see Appendix 
B), the following approximate result up to leading order is obtained, 

1 f°° 2 

Tub = / , (s- N3j)q 2 (s)ds 


N 2 

1 

'W 

1 

N 2 


lN3j 

f(s-N-sj)(^)ds 


“A 79 / (s- Ntj)q 2 (s)ds 


s=N'5 j—>—oo 


s=finite 


l 

'n 2 


, 2 s , e 

( s-N3j )(- 


-#«2 


2^/ns 


r) ds 


(3.26) 


r 

12 ’ 


where we separated the integrand in the first line into three terms for s —> — oo, finite s and 

s —> +oo, since the integration in the limit N — > +oo is from minus infinity to plus infinity. 

2 

Notice that s is bounded from below by AA j and therefore it cannot tend to minus infinity 

2 

faster than N3 j. The third and fourth lines vanish in the large N limit and only the second 
line gives a finite contribution and the result can be obtained by evaluation of the integral in the 
second line. Above approximate calculation produce the correct result as it will be shown in the 
following. 

In the next-to-leading order, the upper bound energy and total free energy of CSM theory 


in each phase can be computed, either by using the above result, Eq. (3.26), or by direct 
computation from Eq. ( 3.25[ ), and the free energy of GWW model, Eq. (3.22), and a new 
obtained result for the expansion of the Tracy-Widom distribution, Eq. ( B.5| ), I IBa-Bu-Dil l. and 
the result is 


(3.27) 


far C < 1 : Fcsm = 



T + N 2 log(l - 


log(c 3 

12 


e - cl j3/2 N 
32nj 3 / 2 N 

e -C2\i\ 3 N 2 N 
|j|l/8ATl/12 ) 


) 


A’ 1 > C + l 
A" 1 < C + l 

























where c\ = 4/3, C 2 = 1/12, C 3 = 2 42e^*)^ ^ and j = X i n this phase, and in the 


other phase, 


(3.28) for C > 1 : Tcsm = 


1 1 

2 _ 3£3 


-C 1 J 3 / 2 N , 


c-!-^ + ^iog(i-S^) A- 1 >2C 


£_ 3 _ 

s 4 


logC 


+ O 2 log(c 3 


e -c 2 \j\ 3 N 2 
|j|i/SAn/i2 


) A " 1 < 2 C 


where j = 2 * , 2C? t t . 

2-3C3(C 2 +C 2 ) 3 

In the large N limit, by expanding the logarithm, the CSM free energy in leading order can 


be computed from above equations, 


(3.29) 


and 


for ( < 1 : Fcsm = 


4 - #r?l A “ (C+ 1)| 


A " 1 > C+l 
A ” 1 < C + l 


(3.30) for C > 1 : TCSM = 


lo §C 
^4 2 



C 2 


^t T |A“ 1 -2C 


2-2C(C2+C" 5 ) 

Thus, in different regions of the moduli space, the upper bound energy is given by 

-^lA" 1 - (C + 1)| 3 C < 1, A” 1 < C + 1 


A " 1 > 2C^ 
A -1 < 2(1 


(3.31) 


•Fuh 


C 2 


2 - 2 C(C 2 +C“ 3 ) 


— IA - 1 — 2 C 2 1 3 £ > 1, A ” 1 < 2( 


After inserting C 2 = 1/12, one can see that the above result is exactly the same as Eq. (3.26 1 . It 


would be interesting to compare the above explicit expressions for the free energy, and possible 
analytic expressions for free energy that would be obtained from Eq. ( A. 1 7[ ) by putting the 
equilibrium eigenvalue density of CSM theory with GWW potential. However, by using Eq. 
( A.17[ ), the general formula of free energy is used in the discussion of level-rank duality, in lHall 
and BTakll . but in comparison to our studies, the analytic formula for free energy as a function 
of A and £, is not obtained and consequently the orders of the phase transitions in the CSM 
theory with any potential such as GWW potential arc not discussed in IITall and l lTak ll. After 
some remarks, based on the obtained results we will discuss the phase structure and the orders 
of phase transitions in CSM theory. 


It can be seen from Eq. (3.31) that, depending on the phase, the upper bound forms at a 
critical coupling A c , which can be determined exactly in each phase, 


(3.32) 


Ac 1 = 


2(l-e)£5 

( 1 - 0(1 + 0 


C> 1 
C < 1 ’ 


and upper bound critical energy, Tfff 1 , infinitesimally close to the domain wall, at A c , is 

32c 2 e 3 £§ / (£§ + £§) 

2c 2 e 3 (l + £) 3 /C 


(3.33) 


-pcrit _ 

** ub 


C> 1 
C< 1 


The free energy of CSM theory, Eqs. ( |3.29 ) and ( |3.30[ ), determines the phase structure of the 
CSM theory, Figure (1). In fact, the explicit formula for the free energy of the upper bound, 
determines the order of the phase transition between the continuous theory (YM theory) and 

13 



























1 


A 



FIGURE 1. Phase structure of Chern-Simons matter theory with GWW poten¬ 
tial. Black line (separating regions II and III) and blue line (separating regions 
I and II), and red curve (separating regions III and IV) are third-order domain 
walls and green line (separating regions I and IV) is a second-order domain 
wall. 


discrete theory (CSM theory). From the upper bound energy, Eq. (3.311, it is easy to find that 
the phase transition at A -1 = ( + 1, between regions (I) and (II), and the phase transition at 
A -1 = 2(2, between regions (III) and (IV) are third order, i.e. the third derivative of the free 
energy, with respect to A or (, jumps at the critical line. By direct computation one can observe 
that a second-order phase transition between regions (I) and (IV) happens at domain wall ( = 1. 
One needs to check the continuity of derivative of free energy, with respect to (, in regions (I) 
and (IV) at C = 1 up to second derivative. However, this domain wall separates also regions 
(II) and (III) and in those regions, this is a third-order domain wall by simple computation 
from YM free energy, Eq. (3.22 1 . Thus, upper bound changes the third-order domain wall to a 
second-order domain wall. 


In summary, the above results in CSM theory with GWW potential determine the phase struc¬ 
ture of the theory and as a new result of our study we have determined the order of phase 
transitions. In this theory, there are four different domain walls separating four different phases. 
I) At the GWW domain wall, the black line, £ = 1, where a third-order phase transition occurs 
and II) at the green line, ( = 1, separating two discrete theories, a second-order phase transi¬ 
tion occurs, and III) two third-order domain walls separating continuous theories and discrete 
theories; the blue line, A -1 = C + lin£<l phase and the red curve, A -1 = 2(2 in ( > 1 
phase. 


3.4. Consistency checks and remarks. Consistency checks can be observed in the following. 
The YM eigenvalue density, p(a), Eq. (A. 19!, is maximum at a = 0. Hence the upper bound 
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condition, p(a) < becomes 


(3.34) 


p(a = 0) 


^-(1 + Ccosa) 


2n 

c 


— + 0 < 2ttX 


Q!— 0 


cos(a/2) Y^~^sin^(a/2) 


= ^< 


Q—0 


2ttA 


C< i 
C> 1 


These conditions are precisely the conditions in F v j, = 0 or TZ(X,() = 1, in regions (II) and 
(III). In this case, the CSM theory is described by the continuum YM theory and thus the 
eigenvalue density of the CSM theory, which in the phase 1Z(X,() = 1 is the eigenvalue density 
of YM theory, should satisfy the upper bound. 

In the following, some explicit values of the upper bound energy in strong and weak coupling 
regime are computed. Taking the strong coupling limit, A —> oo, is only possible in regions (I) 
and (IV) and the values of the upper bound energy at strong coupling are given by 


(3.35) 
and 

(3.36) 


in region (I) : F v }, = 


in region (IV) : F u b = 


16c 2 

C->1 

oo 

C^o’ 

16c 2 

C-» i 

32c 2 

( —y oo 


More generally, in region (IV), for any fixed A, F u b = 32c2 as ( — > oo. In the zero coupling 
limit (free limit), A — > 0, we expect to see the continuous YM theory and thus T u b = 0. This 
limit is only possible and meaningful in region (IV) (because in regions (II) and (III) the CSM 
theory is always in continuous phase) and in this region it implies £ —» oo which leads to 
F u b -A 0, as A -A 0 and ( -A oo, because A^ 1 — 2(2 —> 0, as e — > 0. Such behavior in zero 
coupling limit is in fact consistent with results of other studies in USh-Yil . llGll and UCh-Mil . 
It has been shown in those studies that the thermal partition function of CSM theory in the 
free limit is governed by continuous YM partition function ( |A.1[ ) and the phase structure of the 
theory is similar to that of YM theory. 

On the line A = in region (II), F u b = 0 and in region (IV), F u b — — ^ 32c2 i (1 — (a) 3 - 

C(C 2 +C 3 ) 

At the point ( = 1, A = F u b = 0 and at ( —»• oo, A = F u b = 32c2. 


4. Interpretations and discussions 


There are some possible explanations as well as interpretations of our results. A possible 
explanation is via the understanding of the physical picture of the phase transition. The level 
k is a formal source of discreteness in the system and the real source of discreteness is the 
existence of the upper bound. In fact, since all the parameters are going to infinity, taking 
the continuum limit strictly depends on the relations between the parameters. The plausible 
continuum limit exists under certain conditions. In the naive continuum limit, we expect to 
remove the upper bound, A -» 0, on the eigenvalue density, but at the same time, for larger 
and larger T the eigenvalue density is sharper and peak of the eigenvalue density is higher, so, 
A —> 0 limit should be taken before the limit T —y oo. On the other hand, for finite A in some 
phases, although the eigenvalue density is bounded but the upper bound is not saturated and the 
upper-gap is not formed and therefore the continuum limit of the CSM theory exists. Therefore, 
we can explain our results as a sufficient and necessary conditions for the continuum limit. In 

15 













summary, in order to get the expected continuum limit, k should go to infinity faster than N 
(and T), i.e. the relations k > N + T and k > 2\JNT in two phases. 

In this study we have observed a new second order phase transition between two CSM the¬ 
ories in regions (I) and (IV). It would be interesting to understand possible physical meanings 
and implications of this second order domain wall. 

In the upper-gap phases, the partition function and free energy of the CSM theory, in addition 
to the YM part, contain an extra term originated from the upper bound. The upper gap is a 
consequence of discretization of the eigenvalues of the holonomy matrix and the discreteness is 
itself a direct consequence of the summation over flux sectors in the partition function. Discrete 
fluxes are dual to the monopole operators and discrete flux configurations are effectively the 
same as the configurations generated by monopoles at the centre of S 2 . However, in our case, 
there is no physical monopole or any matter inside S 2 . Having said that, we might interpret the 
formation of the upper bound configuration and its energy T,,}, as the formation of the monopole 
in the CSM theory with the same energy. Possible physical implications of our results in the 
context of monopoles, e.g. monopole condensation and confinement is highly interesting and 
remain for future studies. 

As it is discussed in IfTal and ITakl and references therein, the level-rank duality in CSM 
theory leads to conservation of the partition function under the transformation N —> k — N and 
k —> k or equivalently, A —> 1 — A and ( —> (( A/(1 — A)). At A = A c = 1/2, the level-rank 
duality is valid. But in the general case, the level-rank duality does not seem to be valid. At 
A = A c , the CSM free energy in each phase satisfies the level-rank duality but further studies 
are needed to find other solutions, A*, satisfying, 


(4.1) 


Fcsm^, C) 


= Fcsm{ 1 - A, 


CA , 

1 - A' 


A similar third-order phase transition, between discrete and continuous two-dimensional 
gauge theories, occurs in the Douglas-Kazakov model, [Do-Ka il. Similar phase transitions hap¬ 
pen also in 2d q —deformed YM theory, lArl . IlCal and MJa-Mal . It would be interesting to 
explore possible relations between our study and these works. 

In this paper, the CSM theory is considered with the GWW potential, however, the general 
results in section (3.1) allow to study the asymptotic limit of the ratio of the partition functions 
of CSM and YM theories for any arbitrary potential. One would expect to get similar phase 
structure of upper-gap and no-upper-gap phases for an arbitrary nontrivial potential. We con¬ 
jecture that the orders of phase transitions between different phases of CSM theory with any 
potential in Fig. (1) is universal second and third orders. However, the actual computations of 
the free energy and the order of phase transitions remain for future studies. 

Possible generalizations of the obtained results in this paper, e.g. the third-order phase tran¬ 
sitions, to CSM theory with other gauge groups G = 0{N) and Sp(N), remain for future 
studies, UCh-Joll . 

Another direction for further study is possible applications and implications of our results, 
especially the new phase structure of the model, in AdS/CFT correspondence with Vasiliev 
higher spin gravity. 


Acknowledgements. Thanks to N. Jokela and T. Takimi for useful discussions and reading 
the manuscript. 
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A. Random matrix theory, Toeplitz determinant and Riemann-Hilbert 

PROBLEM IN GAUGE THEORY 

In this appendix, two-dimensional Yang-Mills theory and three-dimensional Chern-Simons 
matter theory, and their relation to matrix models and Toeplitz determinants are reviewed. The 
examples of YM and CSM theories with the Gross-Witten-Wadia potential are considered and 
explicit formulas for the partition function and free energy are reviewed. Furthermore, the phase 
structure of the models and the eigenvalue distribution are described. Finally, the GWW model 
as the solution of the Riemann-Hilbert problem is briefly stated. 

A. 1. YM theory, CSM theory, matrix models and Toeplitz determinant. 


YM theory and CSM theory as matrix models. The interesting and rich phase structure of 
2d YM theory and its matrix model representation with some special potentials have been 
studied in both mathematics and physics literature, HGr-Wil . IWal and HBa-De-Jol . The two- 
dimensional YM theory is defined via the following matrix integral representation of the parti¬ 
tion function, 


Z Y m = / [DU}e 


—Vym — 


(A.l) 


iV COO / 

n / da -j n ( 2 sin ( 

N COO 

n / da ^ n e ' 

• 1 J — OO 7 


^t a l~ a P\\ e -V YM 


im _ „ict p 12 — Vym 


where Vym is the potential of the YM theory. This is a partition function of coulomb gas of unit 
charges on the unit circle with an attractive potential Vym and a logarithmic repulsive potential, 


log | 


ai 


— e 


I, from the Vandermonde determinant. 


Similar to the above result, the path integral of the CSM partition function Eq. (2.8 1 can 
be evaluated by using the techniques from IIBl-Th 1 for pure CS and its modification to include 
the effect of potential, IfTSl . For the study become self contained, we summarize the steps of 
such computation and sketch the derivation of each step but for further details and implications 
we refer to the mentioned references. As we observed in section (2), after integrating out the 
massive fields in the path integral \2A\ , we are left with a path integral ( |2.8[ ) over massless fields 
and the holonomy field. In the following, we integrate on the remaining massless fields and we 
reduce the path integral to a matrix integral over the holonomy fields. 

The first step towards this goal is fixing the gauge invariance of the theory. Let us write A as 
A \ : 2 along the S 2 direction and A 3 along the .S ' 1 direction. Then, there are three gauge fixing: 
i) temporal gauge, d^U = 0 for U = which is the closest gauge to A 3 = 0 and leaves 
a 2d gauge invariance. For further constraining this gauge freedom, we impose ii) diagonalize 
U(x) simultaneously for all x £ S' 2 . This constraint Abelianize the 2d gauge invariance to a 
2d U(1) N gauge invariance. We further constrain the residual 2d Abelian gauge invariance by 
imposing iii) a Coulomb gauge on the time independent diagonal elements of A 12 . 

As we mentioned, in order to evaluate the path integral in ( |2.8| ), we want to integrate over 
the remaining fields Ai^ and reduce the path integral to a matrix integral over U. To do so, we 
impose the above gauge fixing and after some straightforward computations the path integral 
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( |2.8| ) reduces to 
(A.2) 

Zcsm = [ [VA 3 ][VA 2 ][VA 1 }A F pe- i £f d 3 xTr ( A 3 ( dlA 2 - a 2 Al '> +D 3 AlA 2 ) - T J' d 2 x ^ v ( u< 'X'>\ 


where App is the Fadeev-Popov determinant of gauge fixing conditions and D 3 = d 3 + [A3, ]. 

Now by using the Gaussian integration and Lagrange multiplier method it is straightforward 
to integrate over A\^ in the following manner. Notice that, because of the gauge fixing, the first 
term in the exponential receives contributions only from time independent diagonal elements of 
A 1/2 and these elements do not contribute to the second term in the exponential. Then, obviously 
the path integral over off-diagonal elements of A\^ in the second term is of the Gaussian form 
and yields , 1 =, where Dety(D 3 ) is a determinant of the operator D 3 acting in wedge 

y/Det v (D 3) 

product over vector functions. 

By imposing the gauge condition (iii), the path integrals over time-independent diagonal 
elements of A] 2 in the first term of the exponential, evaluated to a delta function that fixes U(x) 
to be constant on S 2 . Then, because of the delta function, the path integral over the eigenvalues 
of the holonomy matrix is reduced to an ordinary integral over the constant eigenvalues of 
holonomy matrix. Lets m/s be constant units of flux for the j— th U( I) factor, then after 
summing over all flux sectors, the full path integral for holonomy matrix reduces to, 


(A.3) 



gi kmjQj 


As we will see, the above contribution from the [7(1 ) A flux sectors leads to the discretization 
of the holonomy eigenvalues. 

It is heuristically argued in IfTSTl and carefully evaluated in llBl-Thl that the Fadeev-Popov 
determinant of the gauge conditions is A pp = Dets(D 3 ), a determinant of the operator D 3 
acting on scalar functions on S 2 and also \J Defy (I)/) = Det' s (D 3 ), which is a determinant 
of the operator D 3 acting on scalar functions on S 2 excluding the zero mode of the scalar field. 
Therefore, the Fadeev-Popov determinant mostly cancels by the path integral of the off-diagonal 
elements and the ratio of Fadeev-Popov determinant and determinant from the path integral over 
off-diagonal modes, which is independent of the flux sector, becomes 


(A.4) 


A F p 

\J dety(L>3) 


detg(£> 3 ) 
det s(D 3 ) 
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Putting above pieces of calculation together, partition function of CSM theory, Eq. (2.8 1 , is 
given by a discrete unitary matrix integral. 


ZcSM = j [DA] e (.iSc S - T f d 2 x^gv(U(x))) 

,cq(rt) - oipilt) N 2 


(A.5) 


iV POO / 

n/ da j\Yi ( 2sin ( 

j=l J -°° Kp^ 

N POO / 

11/ da j []^[ ( 2 sin( 

j=l J -°° Kp^ 

N oo / 

n e o( 2 -( 2 




e ikm J a J 


m, =—oc 


i =1 rii=—o o Z<i 


cq(n) - a p (rt). 


. ,ai(rt) - ai(~rt). 


e ^ 8(kaj — 2irn) 

nEZ 


where from the second line to the third line we used Ylm =-oo e lkm i a .i = 8(kaj — ‘Inn), 

and the delta function puts restriction on eigenvalues, ay(TZ’) = l7V ^ :l , such that no two n, arc 
allowed to be equal, which in return transform the integral over the eigenvalues to a discrete 
sum in the last line. 

A. 1.1. Toeplitz determinant, matrix model and Brownian motion. The relation between matrix 
integral and Toeplitz determinant is briefly introduced in the introduction. In this part, in order 


to elaborate more on Eq. (A.9l, the continuous and discrete Toeplitz determinants are defined 
and their relations to continuous and discrete matrix models arc stated. 

A continuous Toeplitz determinant with a weight function / is defined by 

JV-1 


(A.6) 


Djv(/) = det 


~-j+i 


r l*l=i 


/(*) 


dz 

2i:iz 


J j,l=0 


where z = e ia . Similarly, the discrete Toeplitz determinant is defined in a finite domain d with 
the size |d| as follows 


(A.7) 


D# (/,<*) = det 




-3+1 


\d\ 


m 


z£d 


N -1 


3,1 =0 


The discrete and continuous Toeplitz determinants are related to the discrete and continuous 
matrix integrals via a relation called, Heine-Szego identity, 

N r 2n 


d "(/> = 


(A.8) 


D!?'(/, d) = ' 


i<p 
N 

jviidpv e n/fe) n 

(zi,...,z N )£d N j=l 1 <j<l<N 


Zj - Zi\ 


Therefore, the partition functions of YM theory and CSM theory are equal to the determinant 
ofNxN continuous Toeplitz matrix (continuous Toeplitz determinant), Dy(/), and discrete 
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Toeplitz determinants, D^(/, d\), on domain d\ = {z E C| z k = 1} with size |di| = k, 
respectively, 

(A.9) Z YM = Djv(/), Z CSM = Djy(/, di), 

where / = e~ VyM is the generating (weight) function. 

In IBa-Lil discrete Toeplitz determinants are used to study non-intersecting Brownian mo¬ 
tion. The main result of that study is the following result, Theorem (1.1) in IBa-Lil . Let d be 
a finite discrete subset of a unit circle S 1 with size \d\ and let Q be a neighborhood of S ' 1 and 
f(z ) be an analytic positive weight function with continuous measure, f(z) , then, 

„ D$(/,d) 

<A ' 10) “DM7T =det(1 + A> ’ 


where K is the integral operator with kernel K(z , w) of the Fredholm determinant. And as we 
briefly pointed out, the kernel is given by a factor of continuous kernel as 
(A. 11) 


K(z,W) = Kcont{z, w)y/v(z)v(w)f(z)f(w), 


Kcont^Zi 


nPn{z)p * n ( w ) -P*n( z )pn(w) 
z 1 -z-'w 


and pn(z) is an orthogonal polynomial with respect to weight function f(z), and p* N (z) := 
z n pn{z~ 1 ), and v(z) is a function representing discreteness. This function is defined by 


(A. 12) 


v(z) 


Z7'(z) 

fc7(z) 

nrf W _ | 

k'y(z) 


z G Sf n 
z G Sl M ’ 


where 7 (z) is an analytic function on 12 such that it vanishes exactly on d and all zeros are 
simple, and Sf n and Sg Ut are positively-oriented circles of radius 1 — e and 1 + e, for small e > 0. 
For proof of the above result in a slightly different notation, see llBa-Lill . 

Toeplitz determinant appears naturally in the context of Brownian motion as probability. In 
order to state some other useful and relevant results from HBa-Lil which we are going to use 
them in gauge theory, let us briefly introduce a minimal necessary notation in random walk 
theory. Let Xi(t), i = 1..... A r be independent standard Brownian motion with conditions, 
Xq( t) < X\{t) < ... < .Xjv-i (t) for all t G [0, T] and 2Q(0) = Xi(T) and define the width 
by Wn{T) = sup te [ 0 7 ’] (An- 1 (t) — X\ (t)). Then, the conditional probability for the width 
Wn(T) of non-intersecting continuous time symmetric simple random walks is obtained in 
(Proposition 4.1), IIBa-Lill . 


(A.13) 


P (W n (T) < M) 


lim 

N,T,k^-oo 



^N^ifoww, d s ) ds 

D n(Igww ) 27tls’ 


where d s = {z E C| z l = s}, and the asymptotic limit is obtained in (Theorem 4.1), liBa-Liil . 


(A. 14) 
where 

(A.15) 


lim P ( 

7V,T->oo 


W n (T) — p(N, T) 
a (N , T) 


<x) = F(x), 


F := 


N + T N >T 
2 y/NT N <T 


a := 


1 1 

3T3 




and F(x) is the Tracy-Widom distribution, discussed in Appendix B. 
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A. 1.2. Free energy and eigenvalue distribution function. In the YM theory, in order to find the 
equilibrium measure one has to look for the measure on the circle which is non-negative and 
satisfies an energy minimization problem. 

The distribution of the eigenvalues on the unit circle is determined by solving the follow¬ 
ing variational problem, HBa-De-Joll . For a given measure f(z) yyf, the equilibrium measure 
(eigenvalue density) p(z) is obtained as a unique minimizer of the free energy in following 
variational problem, 

(A.16) E ' = infjj 7 ^ ((p)) : pis a probability measure on the unit circle S' 1 }, 

where 

(A.17) E x (p) = / I log \z — w\~ 1 dp(z)dp(w) + / V(z)dp(z), 

J Js 1 xS 1 Js 1 

is the free energy. 


Gross-Witten-Wadia model. Consider a specific potential Vym = — (T/2)Tr(U+U^) where U 
is a N x N holonomy matrix with eigenvalues o:,’s. This defines an interesting two-dimensional 
YM theory, called Gross-Witten-Wadia model, IlGr-Will . IIWall , which is a unitary random one- 
matrix model. The partition function reads, 


~(GWW) 
^YM 


(A. 18) 


f [VU]e^ u+u ^ 




T 

e 2 


( z j+ z j x ) 



_ e ia p 12 e T J2f =1 cos aj _ 


The above partition function is equivalently given by the Toeplitz determinant, D n(/) with 

j _ gT cos a 

By solving the variational problem for the GWW model, V(z) = ^(z + z -1 ), the eigenvalue 
density is obtained, 


(A. 19) 


p(a) 


A;(l + Ccos a) (<1 

£cos(a/2)^ - sin 2 (a/2) (>1 


As it will be explained in below, this model possesses a third-order phase transition and the 
phase transition happens at £ = 1 , and the eigenvalue density p(a) has different features in 
different phases. In ( < 1 phase, p(a) is supported on whole circle and this phase is called no¬ 
gap phase. In £ > 1 phase, p{a) is supported on a subset of circle, [— a c , a c ] where a c is given 
by the condition sin 2 (a c /2) = 1/C and this phase is called lower-gap phase. Free energy of the 
GWW model is defined by = limjv^oo 777 log Z^W W \e ( ’ N ^-' : > cosaj ), and can 

be computed by putting the eigenvalue density (A. 19) in ( |A.17[ ) and with a change of variable 
T = C-A, it has been computed in I Gr-Wil . 


(A.20) 


t (gww) 
j "ym 


(0 
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A third-order phase transition, i.e. discontinuity in d 3 F/d 3 (, occurs at £ = 1. 


CSM theory with GWW potential. In the following, the qualitative and quantitative results for 
the phase structure of CSM theory with GWW potential are reviewed. The explicit results 
are obtained from the saddle point analysis of the variational problem for the discrete GWW 
model, in low-temperature H'Tal and high-temperature UTakl . For alternative derivations of the 
eigenvalue measures for models of this type, see HJo-Ja-Kel . 

At low temperatures, in the regions (A < 1/2, £ < 1) and (A > 1/2, £ < 1/A — 1), the 
repulsive force (factor (II) in discussion in section (3.2)) is stronger and eigenvalue density 
has support everywhere on the unit circle. The eigenvalue density is given by pi a) = (1 + 
C cos ol)/2-k and the system is in no-gap phase. As the temperature increases, the attractive force 
(factor (I)) plays more effective role and depending on A, the lower-gap or the upper-gap forms 
first. For small A, A < 1/2, in the region 1 < C < 1/4A 2 , the upper bound is large and first 
the lower-gap phase forms with eigenvalue density, p(a) = C, cos (a:/2) /tt \/l/C — sin 2 (a/2) 
supported on [— a c , a c ] and a c is given by the condition sin 2 (a c /2) = 1/C- Then, in the region 
C > 1/4A 2 the two-gap phase with one lower-gap and one upper-gap forms. The lower-gap and 
the upper-gap are located on the arcs [ e lb , e~ lb ] and [e ia , e~ ia ] on the unit circle, respectively, 
where a and b are determined by 


= c, 


1 f a . cos a 

(A-21) 1—\ / da I 7= 

47rA •'-» \J sin 2 | — sin 2 ^ys 


l 2 f\/sin 2 ! -sin 2 f 


1 + ^(cos ° + cos b). 


Between the arcs, in the complement of two gaps, the eigenvalue density is determined by 


p{a) = 


(A.22) 


| sin a | 
47t 2 A 

/ dO 


f Q OL n CL . . r\ b r\ Oi . 

.sin 2 _ sm ^ 2 Sin 2 S11 ^ 2 ’ 


(cos 6 — cos a) J sin 2 % — sin 2 | J sin 2 | — sin 2 ^ 


For large A, A > 1/2, in the region 1/A — 1 < C < 1/(4A(1 — A)), the upper bound is small 
and upper-gap phase forms first, with eigenvalue density, 


(A.23) 


p(a) = 


1—A 
AC 


9 ol 
— COS z 2 


cos 2 § < 


cos 2 f > 


AC 

1—A 

AC 


then, in the region C > 1/(4A(1 — A)), the two-gap phase forms. 

In summary, the phase structure of the CSM theory with GWW potential is 
(A.24) 


for A < A c : < 


no-gap phase 

C< 1 


no-gap phase 

c<i-i 

lower-gap phase 

1 < C < 4^2 ) 

for A > A c : < 

upper-gap phase 

A-!<C< 

two-gap phase 

£ > D? 


two-gap phase 

£ > 4A(1—A) 


1 

4A(1—A) > 


where A c = /. In principle, the above results for eigenvalue density can be put in Eq. (A. 17 1 
to give the free energy of GWW model. The general form of free energy, Eq. (|A. 1 7|), has been 
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discussed and used for CSM theory with different potentials including the GWW potential in 
the discussion of level-rank duality, llTakl . However, our new explicit results for free energy 
of the CSM theory with GWW potential as an analytic function of A and £ is obtained from a 
completely different method. 


A.2. Riemann-Hilbert problem in Yang-Mills theory. In this part, the YM theory as a con¬ 
tinuous matrix model is formulated via a well-defined mathematical problem, namely the Riemann- 
Hilbert problem (RHP), IIBa-De-Joll . In the following, orthogonal polynomials and RHP asso¬ 
ciated to the GWW model are reviewed. Lets zGC, and T E 1R is a parameter of the model, 
then the following polynomial, 

(A.25) p N (z,T) = a N (T)z N + ... 

is called an orthogonal polynomial with respect to weight function, f(e ia ) = exp ({T /2)(e ia + e~ ia )), 
on the unit circle if it satisfies 

/ 7T _ 

PN(e ia )pM(e ia )dfj,f = 5 n,Mi 

-7T 

where N,M > 0 and measure is d/ij = /(e 1Q )|A The well-known connection between the 
orthogonal polynomial and Toeplitz determinant is established in [[SzJ. As a result, the leading 
coefficient in orthogonal polynomial, a/v(T), is expressed in terms of the Toeplitz determinant, 

, A frr ^ Dtv_i(T) 

(All) a N (T) = 

where B N (T) = B N (e Tcosa ). 

It has been observed earlier that the partition function of the YM theory with the GWW 
potential which is given by a unitary matrix integral is expressed in terms of the Toeplitz deter¬ 
minant. Thus using the above connection to the orthogonal polynomials and the Szego strong 
limit, lmpv_>.i D n(T) = 1, one can obtain, 


(A.28) log 4^ } (T) = logD*_i(e Tcosa ) = Y, kg a?(T). 

l=N 


To control the behavior of the partition function in the limit of large N, T, one needs to control 
the behavior of orthogonal polynomial af(T) in the large N. T limit for all l > N. To this 
aim, the steepest descent methods has been used to compute the asymptotic behavior of the 
following RHP, IIDe-Zhll . 

The following RHP uniquely determines af(T) and partition function of the YM theory 
with the GWW potential. Let S be a unit circle with the counter clock-wise orientation and 
Y(z, l + 1, T) is a 2 x 2 matrix-valued function satisfying 


'Y(z, l + 1, T) is analytic in C — £, 

< Y+(z,l + l,T) = Y-(z,l + l,T) 

Y(z, l + 1, T)z-^ +1 ^ 3 =1 + 0(1) asz^oo 
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1 ~kre^ z+z ^ 


(A.29) 










where Y± denotes Y inside/outside of the unit circle, 03 = 


1 0 
0 -1 


and z ( i+1 ) 0 ’ 3 = 


z ~(l+i) Q \ 

l+1 J . It has been proved that the largest coefficient in the orthogonal polynomial 
is given by 21—entry of matrix Y, 

(A.30) a i(T) = — l2t(0; l + 1, T). 

Asymptotic results for a 2 N (T) and 1 I'or A". T —> oo is obtained in HBa-De-Jofl . 

B. Tracy-Widom distribution 

Tracy-Widom distribution Fix), ll Tr-Wil . HBa-De-Jol . is the probability distribution of the 
largest eigenvalue of a Hermitian random matrix and it is defined by 

(B.l) F(x) =e-f*°° (s - x)r}2 ^ ds , 

where q(s ) is the solution of the Painleve II equation, q"(s) = 2 q 3 (s) + sq(s), with boundary 
condition q(s) ~ — Ai(s ) as s -X oo, where Ai(s) is the Airy function. 


(B.2) 


Ai(x) = A 


. 3 1 r oo 6 

e J t +lxs ds = — cos( + xs)ds. 

-oo ^ JO ^ 


The Tracy-Widom distribution is also given by F(x) = det(l + Ka,') where the kernel of Airy 
function is given by 

Ai(z)Ai'{w ) — Ai(w)Ai'(z) 


(B.3) 


K A i(z,w) = 


z — w 


In the asymptotic regime, s —> Too, the following result is obtained in llBa-Bu-Dil . 


(B.4) 


q(s) = < 


4 ! 


-Ai(*) + 0(M-) 


1 

S4 


Too 


-Afl + _i_Z3_ j_ 

2 dr 8s 3 !28 s g ^ 1024s 9 


1020 


+ d(|s| 12 )) s -> —oo 


- 2 S ! 
e 3 

2y/7TS^ 


where in the limit s —>• +oo, Ai"(s) = sAi(s) and 

Moreover, the asymptotic of the Tracy-Widom distribution is studied in BBa-Bu-Dil . and the 
following result is obtained, 


(B.5) F 0) = < 

where is the Riemann zeta function. 


1 - 


4 § 
e-a * 3 


*-( 1 -^ 3 - + 0 ( 0 ) 

32?rx2 24a; 2 

-Ahl 3 


2i e «-«-DO^L(l + 3 +0(W -6 ) ) 

\X\ 8 


as x —}■ oo 
as x —> —oo 
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